in 
o 

o 



Particle Creation by a Moving Boundary with 
Robin Boundary Condition 



B. Mintz*, C. Farina 1 ", P.A. Maia Neto* and R.B. Rodrigues § 



O 
O 

(N 

J>^! Instituto de Fisica - Universidade Federal do Rio de Janeiro 

Caixa Postal 68528 - CEP 21941-972, Rio de Janeiro, Brasil. 



(N 

Abstract 

We consider a massless scalar field in 1+1 dimensions satisfying a 
Robin boundary condition (BC) at a non-relativistic moving boundary. 
s ! , We derive a Bogoliubov transformation between input and output bosonic 

field operators, which allows us to calculate the spectral distribution of 
created particles. The cases of Dirichlet and Neumann BC may be ob- 
tained from our result as limiting cases. These two limits yield the same 
spectrum, which turns out to be an upper bound for the spectra derived 
for Robin BC. We show that the particle emission effect can be consider- 
ably reduced (with respect to the Dirichlet/Neumann case) by selecting a 
particular value for the oscillation frequency of the boundary position. 

1 Introduction 

Moving bodies experience fundamental energy damping £Q [2] and decoherence 
r> | |3] mechanisms due to the scattering of vacuum field fluctuations. The damping 

■ is accompanied by the emission of particles (photons in the case of the elec- 

tromagnetic field) 0|, thus conserving the total energy of the body-plus- field 
system This dynamical (or nonstationary) Casimir effect has been an- 

alyzed for a variety of three-dimensional geometries, including parallel plane 
plates |7|, cylindrical waveguides and rectangular [S], cylindrical ^U] an d 
spherical cavities It also depends on the details of the coupling between 

the field and the body, which can usually be cast in the form of boundary con- 
ditions (BC) for the field. Of particular theoretical relevance is the Robin BC, 
which continuously interpolates the Dirichlet and Neumann BC. For a massless 
scalar field in 1+1 dimensions, it reads 

dcj) . . 1 , . . . 

^(^o) = ^Wi a; o), (1) 
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where xq is the position of the boundary. The positive parameter (3 represents 
a time scale (we take c = 1) associated to the time delay (or phase shift) 
characteristic of reflection at the Robin boundary which can be interpreted 
in terms of a simple mechanical model ^3]. According to eq. (Q, Dirichlet and 
Neumann BC are obtained as the limiting cases j3 — > and f3 — > oo, respectively. 

We consider a semi-infinite slab (extending from — oo to x = Sq(t)) following 
a prescribed nonrelativistic motion, with the Robin boundary at 5q(t). We have 
recently computed the dynamical Casimir force on the slab, which contains 
dissipative as well as dispersive components |12j . In this paper, we will analyze 
in detail the particle creation effect and compute the corresponding spectral 
distribution. 



2 Input-Output Bogoliubov Transformation 



In the instantaneously co-moving frame the massless scalar field satisfies 

Ml -I^/i 

| bound — I bound i 

Neglecting terms of the order of [(5g(t)] 2 , we find, in the laboratory frame, 



dx +Sq{t) di 



4>(t,6q(t)) = -4>(t,6q(t)). 
P 



(2) 



(3) 



We assume that final position coincides with initial one, which is taken at 
x = 0. Hence 

lim Sq(t) = 0. (4) 

t — >±oo 

Jointly with the nonrelativistic approximation, this condition implies that 8q(t) 
is much smaller than the wavelengths A of the created particles. In fact, we 
will show that the frequencies of the particles are bounded by the mechanical 
frequencies ojq : lu = 2ir/\ < luq. Since u>oSq ~ 6<j <C 1, we have Sq(t) <C A. 

Thus, we may analyze eq. © by expanding up to first order in Sq and 
its derivatives. This amounts to calculate the effect of the motion as a small 
perturbation 

<p(t,x)=<f> (t,x) + 6(f>(t,x), (5) 

where the unperturbed field tf>o corresponds to a solution with a static boundary 
at x = 0. The first-order field 8(j> then satisfies the following BC at x = : 



^(t,o)-ity(t,o) = a«(t) 



6q(t)-^-(t,0). (6) 



It is convenient to use the Fourier representation 
$(w,£b) = / dte luJt cj)(t,x). 



The unperturbed field satisfies the Robin BC at x — 0. Its normal mode expan- 
sion for x > is given by 

<&o(to,x) = N(u)[sin(ujx) + u/3 cos(lox)} [Q(oj)a(u>) - 6(-w)a(-w) t ] , (7) 
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with 



N{lu) = 



4tt 

H(l + (3 2 uj 2 ) 



and Q(x) denoting Heaviside step function. The bosonic operators a(u)) and 
a(ujy satisfy the commutation relation 



[a(uj), a(u/) f ] = 2%6(u-u}'). 



(8) 



To solve eq. © for 6&(u>, x) in terms of $o(w, 0) (with x > 0), we use suitably 
defined Green functions, obeying the differential equation 

d 2 



dx 2 



+ J 2 G(oj, x, x') = S(x - x'). 



(9) 



^From Green's theorem, we find 

6$(w,x') = -5$(uj,0)—G(lj,Q,x') + G(u,Q,x')—8$(l<j,0). (10) 
ox ox 

This result is more easily combined with eq. (0) if we select a solution 
Gr(w, 2;, a;') of eq. © satisfying the Robin BC at x = 0. With this Robin Green 
function, we immediately obtain the first-order field from eq. (|10fl in terms of 
the BC satisfied by 5Q(ui, x) as given by the Fourier transform of eq. ©■ Then, 
the complete field is written as 



$(w, x) = $ (w, x) + G R (w, 0, sc) 



with 
0_ 

dx 



5$(w,0) 



<5$(u;,0) 



/3 
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9£o 

9a; 



(w,0) +ww'$o(w,0) 



(11) 



(5Q(w - a/), 

(12) 

where 8Q{uj) is the Fourier transform of 5q{t). 

If we replace Gr in eq. I|ll|) by the retarded Robin Green function, given by 



(13) 



1 - i/3u 

then the zeroth-order field &o(u,x) corresponds to the input field $i a (uj,x), 
with 

4>i a {t,x) — lim <j)(t,x). 

t — >— 00 

On the other hand, when taking the advanced Robin Green function, given by 







(14) 



1 + 0UJ 

$o(w, x) corresponds to the output field <& out (cc>, x) (</> ut(*> x) = lim^oo cj>(t, x)). 
By combining these two possibilities, we find the relation between output and 
input fields: 



1) = Mw, as) + [GS*(w, 0, 53) - G£ v (w, 0, a)] 



(15) 
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The final result is obtained by inserting eqs. (|T2|) (with $ replaced by $j n sin- 
ce we neglect terms of second order), (fH5)l and (tl4ll into the rhs of eq. (fL5|) . 
Further physical insight is gained if we write the input-output relation in terms 
of the annihilation operators a- m , a out and their Hermitian conjugates, by com- 
bining eqs. (0) and (|15|) . The resulting input-output relation has the form of a 
Bogoliubov transformation: 

2iJu f duj' 1 + 3 2 uu)' r— - 



y/l + fi 2 Lu 2 J 2tt ^1+^72 
x [9{uj')a in {uj r ) - e(-u')a ia (-u>'y]6Q(oj - J). 

Since the output annihilation operator is contaminated by the input creation 
operator, the input vacuum state |0i n ) is not a vacuum state with respect to 
the output operators. In the next section, we compute the resulting particle 
creation effect. 



3 Frequency spectrum 

The number of particles created with frequencies between to and to + diu (lo > 0) 
is 

——{u)du> = (0 in |a out (w) t a ou t(w)|0 in ) — (17) 
dio 2n 

The spectrum is obtained by inserting eq. (|16fl into i|17fl : 

dN, , 2lo f°° duj' (3 2 iolo'] 2 .... , . 

duj y ' ir(l + /3 2 lu 2 ) J 2tt 1 + /3 2 uj' 2 L v n y ' 

To single out the effect of a given Fourier component of the motion, we take 

Sq(t) = Sqo cos(u;o£)e - '*'/ T 

with loqT 3> 1. In this case, 5Q(lo) corresponds to two very narrow peaks around 
lo = ±ujq, so that we may take the approximation 

\5Q{uo)\ 2 « I Sq 2 T [S(co - ujo) + 5{uj + lo )} . (19) 

Inserting this equation into JTHJl, we find 

dN > Sq 2 T [1 - PMcoq - lo)} 2 

Note that the spectrum vanishes for lo > loq : no particle is created with fre- 
quency larger than the mechanical frequency. Field modes at higher frequencies 
are not excited by the motion, which is slow in the time scale corresponding 
to such frequencies (quasi-static regime) . This important property confirms the 
consistency of our perturbation approach, with its expansion in Sq/X. 

A second important general property of the spectrum given by eq. (|2U)l is the 
symmetry around u = loq/2 : the spectrum is invariant under the replacement 
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Figure 1: Spectral distribution of the emitted particles dN/dw. For the horizontal scale, we 
divide the frequencies by the mechanical frequency wq. Dashed line: (3 = (Dirichlet case), 
solid line: /3luo = 1-7, dotted line fluig = 5. 



lo — > ujq — lo. This is a signature that the particles are created in pairs, with 
frequencies such that their sum equals loq. Hence, for each particle created at 
frequency lo, there is a 'twin' particle created at frequency loq — lo. 

Since Robin BC interpolate Dirichlet and Neumann ones, we may derive the 
spectra for these two cases by taking appropriate limits of eq. (|2L)[1 . For Dirichlet 
BC, we find 



dN 



dN . . (5q a ) 2 T . 
= hm—{u) = — u(u>o - w) 0(w o - w), (21) 

, D \ P^O dlO ill 



in agreement with Ref. For the Neumann BC (/3 — > oo), we find the same 
spectrum, confirming the equivalence between Dirichlet and Neumann in the 
context of the dynamical Casimir effect in 1+1 dimensions |14) . 

For intermediate values of j3, the spectrum is always smaller than in the 
Dirichlet case for all values of w. In fact, we may write the result of eq. (|20|l as 

where the reduction factor rj < 1 is a function of [3ujq and u)/u)q. The reduction 
may be more severe near to = u)q/2, which is the spectrum maximum in the 
Dirichlet case, for some values of (3u>$. Hence, the Robin spectrum may develop 
global maxima near lo = and to = ujq, as in the example shown in Fig. 1 
(solid line), with /3loq — 1.7. In this figure, we also plot the Dirichlet/Neumann 
spectrum (dashed line) and the Robin spectrum for f3cvo — 5 (dotted line). As 
discussed above, all curves are symmetric with respect to u) = loq/2. 

The areas below the curves shown in Fig. 1 correspond to the total number of 
created particles. The figure already indicates that this number, to be discussed 
in the next section, may be considerably reduced (with respect to the Dirichlet 
case) for intermediate values of (3. 
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4 Particle Creation Rate 



The total number of created particles is given by 

N = n ^L (uj) du;^ 6 -^ ^ F(Pu; ), (23) 

Jq OUJ Z7T 

with 

p(n m + e + 12 arctan(g)] - 6(2 + g 2 ) ln(l + £ 2 ) 

m ~ wiMe) ( } 

As expected for an open geometry (with a continuum of field modes), N is 
proportional to the time T, so that the particle creation rate R = N/T is the 
physically meaningful quantity. For the Dirichlet case, we take F(£ — ► 0) = 1/6, 
and then 

12tt 

For the Neumann case, we find the same result for the creation rate, since 
the spectrum is the same. Note that the rate increases with too according to 
eq. I|25|) and vanishes (as required) in the static limit luq — 0. This could have 
been anticipated since the particle creation is an effect of changing the BC 
nonadiabatically. However, for Robin BC the rate is not a monotonic function 
of luq. In Fig. 2, we plot the rate R as a function of (3ujo (for a fixed (3). R 
decreases as uj varies from 1.3//3 to the local minimum at 2.1//3. 




Figure 2: Total particle creation rate as a function of mechanical frequency (in units of 1//3). 
Inset: ratio between creation rates for Robin and Dirichlet BC. 



In the inset of Fig. 2, we plot the ratio R/R^d) = 6F(/3ujo) < 1 as a function 
of /3u>o. This ratio represents the reduction of the Dirichlet creation rate for 
a finite p. It only depends on the the dimensionless variable /3uiq, and goes 
asymptotically to one for (3ujq ^> 1, since the Neumann BC yields the same rate 
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as the Dirichlet case. The reduction is maximum at (3uq =2.2. At this point, 
the creation rate is reduced to 1.3% of the Dirichlet value . 

We may also calculate the radiated energy from these results. Thanks to the 
symmetry of the spectrum around lo = ujq/2, we have 

E= n^^Kuidcu^^N. (26) 
Jo " w 2 

Combining with eq. I|23[) . we find 

E = 8q£Ttiu$F(pu )/{4n). (27) 

This expression can be directly compared with the result for the Casimir 
force we have recently reported |12j . The force is written (in the Fourier domain) 
as !F(w) = x(w)5Q(oj) and its work on the slab is given in terms of the imaginary 
part of the susceptibility function x( to ') : 

1 f°° 

W = / rfwwlmx(w) \SQ(lo)\ 2 . (28) 

k Jo 

For the quasi-sinusoidal motion considered here, we find, using eq. (|19(l . 

W = ~Sq 2 Tu Im X (u )/2. (29) 

The result for Im^(wo) derived in Rcf. 12 can be cast in the form 2 Imx(wo) = 
?iujqF((3ujo)/{2tt). Then, the comparison of eqs. and yield E = -W, so 
that the total radiated energy coincides with the negative of the work done on 
the slab by the Casimir force, as expected from energy conservation. 



5 Conclusion 

Dirichlet ((3 — > 0) and Neumann BC ((3 — ► oo) yield the same result for the 
spectrum of created particles. With the Robin BC, we are able to interpo- 
late continuously between these two cases. For intermediate values of j3, the 
spectrum is always smaller than the Dirichlet/Neumann case, for all values of 
frequency. 

In the range 1.2 < (3cuo < 2.4 the spectrum develops lateral peaks higher 
than the value at u> = u>o/2. This is also approximately the range in which the 
total creation rate (surprisingly) decreases with u>q. This rate is reduced by up 
to 1.3% of the Dirichlet/Neumann value, if the mechanical frequency is selected 
at wo = 2.2/ p. In other words, the coupling with the vacuum field state is 
considerably reduced if the slab oscillates at a frequency close to this value. 

When considering the electromagnetic field and a plane mirror moving along 
its normal direction, the BC in the ideal case of perfect reflectors may be decom- 
posed into Dirichlet and Neumann BC for each orthogonal polarization |15|. In 

1 When plotting the creation rate itself, the effect seems to be less impressive because the 
Dirichlet rate increases with loq. 

2 ln Ref. 1121 . a narrow plate is considered, rather than a slab. The two sides of the plate 
provide identical (and independent) contributions to the force. Hence, to compare with the 
present situation, we divide the result of II 2| by two. 
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3+1 dimensions, the effect with Neumann BC is considerably larger than with 
Dirichlct BC, and it would be interesting to investigate the continuous transition 
between these two limiting cases. 

Reflection by real metallic plates involve non-trivial frequency-dependent 
phase factors as in the case of Robin BC. The results of the present paper indi- 
cate that finite conductivity might yield a significant reduction of the magnitude 
of the dynamical Casimir effect. 

This work was supported by the Brazilian agencies CNPq and FAPERJ. 
PAMN acknowledges Instituto do Milenio de Informagao Quantica for partial 
financial support. 
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